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ABSTRACT 

We study the influence of X-rays on the wind structure of selected O stars. For this purpose we 
use our non-local thermodynamic equilibrium (NLTE) wind code with inclusion of additional 
artificial source of X-rays, assumed to originate in the wind shocks. 

We show that the influence of shock X-ray emission on wind mass-loss rate is relatively 
small. Wind terminal velocity may be slightly influenced by the presence of strong X-ray 
sources, especially for stars cooler than T e g < 35 000 K. 

We discuss the origin of the Lx/L ~ 10~ 7 relation. For stars with thick wind this relation 
can be explained assuming that the cooling time depends on wind density. Stars with optically 
thin winds exhibiting the "weak wind problem" display enhanced X-ray emission which may 
be connected with large shock cooling length. We propose that this effect can explain the 
"weak wind problem". 

Inclusion of X-rays leads to a better agreement of the model ionization structure with 
observations. However, we do not found any significant influence of X-rays on P V ionization 
fraction implying that the presence of X-rays cannot explain the P V problem. 

We study the implications of modified ionization equilibrium due to shock emission on 
the line transfer in the X-ray region. We conclude that the X-ray line profiles of helium-like 
ions may be affected by the line absorption within the cool wind. 
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1 INTRODUCTION 

The evolution of hot stars seems to be strongly influenced by the 
loss of a substantial part of their mass via winds. Consequently, a 
proper knowledge of amount of material expelled from the stellar 
surface per unit of time (mass loss rate) is necessary for the stellar 
evolution calculations. 

There exist some serious problems with mass-loss rate predic- 
tions for hot stars. The situation is not satisfactory even in the O 
star domain, where several different theoretical models are avail- 
able (e.g., Vink et al. 2000, Pauldrach et al. 2001, Puis at al. 
2003, Krticka & Kubat 2004, hereafter Paper I). A good agree- 
ment between these theoretical models and observations may be 
just illusory effect of the possible manifestation of local wind in- 
homegeneites (clumping) in the observed spectra (e.g., Bouret et 
al. 2003, Martins et al. 2005, Puis at al. 2006b). The problem is 
that the spectrum of a clumped wind with small mass-loss rate may 
mimic the spectrum of a smooth wind with large mass-loss rate. 
Consequently, the mass-loss rates derived from observations with 
account of clumping may be significantly lower than those derived 
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without taking into account the influence of clumping. Lower mass- 
loss rates are indicated also by weak P V wind lines (Fullerton et al. 
2006), and by the study of X-ray line profiles (Owocki & Cohen 
2006). If the mass-loss rates derived from observations are really 
much lower due to the influence of clumping on the observed spec- 
tra, then there exists a significant discrepancy between theory and 
observation (Puis at al. 2006a). 

Moreover, some O stars with low luminosities exhibit the so- 
called "weak wind problem" (Bouret et al. 2003; Martins et al. 
2004; 2005; Krticka 2006). Theoretically predicted mass-loss rates 
of these stars are much higher than the mass-loss rates inferred from 
observations (even without taking into account the clumping). The 
origin of this discrepancy is unclear. 

One of the reasons that may cause these mentioned differences 
between observation and theory of hot star winds is the neglect of 
X-ray radiation. The X-ray radiation of hot stars is likely generated 
by the wind shocks that exist in the wind either due to the instabil- 
ity caused by the radiative driving (Owocki et al. 1988, Feldmeier 
et al. 1997b) or due to the collisions of wind streams channelled by 
the magnetic field (Babel & Montmerle 1997, ud-Doula & Owocki 
2002). The X-rays may be generated also due to decoupling of wind 
components accompanied by frictional heating (Kubat et al. 2004). 
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In the case of binaries the X-rays may also originate in wind colli- 
sions (Prilutskii & Usov 1976; Antokhin et al. 2004). 

X-ray radiation influences above all trace ionization states in 
most of hot O stars (MacFarlane et al. 1994) and the inclusion of X- 
ray sources is necessary especially for the correct prediction of hot 
star spectra (Pauldrach et al. 2001). On the other hand, for cooler 
stars the wind ionization balance may be dominated by the influ- 
ence of X-ray radiation (MacFarlane et al. 1994). 

Recently Fullerton et al. (2006) showed that the mass-loss rate 
determined from the P V lines is much lower than that derived from 
the Ha line or radio emission. This implies also that hot star wind 
mass-loss rates determined from P V lines are significantly lower 
than the theoretical ones. However, this conclusion is sensitive to 
the ionisation state of phosphorus. If P V is not a dominant ion at 
any effective temperature studied by Fullerton et al. (2006, e.g., due 
to the presence of X-rays), then the discrepancy between theory 
and observation of OB star winds may not be so significant. Hence, 
the inclusion of X-ray ionization may also be important from this 
perspective. 

Thus, we present here a study of the influence of the X-ray 
radiation on the wind properties of O stars. 



2 NLTE WIND MODELS 

The process of calculation of models used in this paper was de- 
scribed in Paper I, thus here we only summarise its basic features 
and describe improvements. 



2.1 Basic model description 

Our models assume spherically symmetric stationary stellar wind. 

Excitation and ionization state of elements important for the 
radiative driving and for the correct calculation of the radiative field 
is derived from the statistical equilibrium (NLTE) equations. 

The radiative transfer equation is artificially split (as in the Pa- 
per I) into two parts, namely the radiative transfer in continuum and 
the radiative transfer in lines. The solution of the radiative transfer 
in continuum is based on the Feautrier method in the spherical co- 
ordinates (Mihalas & Hummer 1974 or Kubat 1993) with inclusion 
of all free-free and bound-free transitions of model ions, however 
neglecting line transitions. The radiative transfer in lines is solved 
in the Sobolev approximation (e.g., Castor 1974) neglecting con- 
tinuum opacity and line overlaps. 

The radiative force (calculated in the Sobolev approximation 
after (I.25) 1 using data extracted in 2002 from the VALD database 
(Piskunov et al. 1995, Kupka et al. 1999)) and the radiative cool- 
ing/heating term (we use the electron thermal balance method, 
Kubat et al., 1999) are calculated using occupation numbers de- 
rived from statistical equilibrium equations. Finally, the continuity 
equation, equation of motion, and energy equation are solved and 
consistent wind velocity, density, and temperature structure are ob- 
tained using iteration procedure. For our calculations we use As- 
plund et al. (2005) solar abundance determinations. 

The lowest wavelength considered in the models is 7.5 A and 
the largest X-ray wavelength is defined as lOOA. 



2.2 Inclusion of shock X-ray radiation 

The consistent inclusion of X-ray generation into NLTE wind mod- 
els (using hydrodynamical simulations that are able to predict the 
shock properties, Feldmeier et al. 1997b) is likely beyond the possi- 
bilities of present computers. To make the problem more tractable, 
and since we are interested in the effect of already generated X-rays 
on the wind, not in the process of X-ray generation itself, we use a 
simpler approach after Pauldrach et al. (1994), i.e., we include X- 
rays into our stationary models in an artificial way using two free 
parameters introduced below. 

We assume that a part of wind material is heated to a very high 
temperature Tx due to the shock. The shock temperature is given 
by the Rankine-Hugoniot shock condition 
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where k is the Boltzmann constant, itih is the hydrogen mass, a» is 
the sound speed calculated assuming a completely ionised hydro- 
gen plasma, 
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where T H is the hydrogen temperature taken from stationary wind 
models, and ux is the upstream shock velocity, 

Ux=U K \V r , (3) 

where v r is the radial velocity in the stellar rest frame and it re i 
is a dimensionless free parameter influencing the hardness of X- 
rays. Instead of the linear dependence of ux on the wind velocity 
it would be possible to employ different assumptions, e.g., con- 
stant ux- However, the observed X-rays originate at larger radii 
where the wind velocity approaches the terminal one, consequently 
such dependence would not result in a very different results. More- 
over, in the case of constant ux it would be necessary to introduce 
additional free parameter corresponding to the radius at which X- 
rays start to be emitted. Without it we would obtain unrealistically 
strong X-ray source at the wind base as the wind density is the 
largest there. 

We add the X-ray emissivity 



1xM = /x(W A„(Tx)/(4jt) ; 



(4) 



into the emission coefficient, where fx is the second free param- 
eter determining the amount of X-rays (also called the filling fac- 
tor), A„(Tx) is calculated using the Raymond-Smith X-ray spectral 
code (Raymond & Smith 1977, Raymond 1988), and £ relates the 
wind density p and the electron number density 

m H 1 + 4y 

where y is the helium number density relative to the hydrogen one 
(here we assume a fully ionised gas). 

We used different values of fx in our calculations. The free 
parameter M re i = 0.3 in all presented calculations. This roughly 
gives the same X-ray temperature compared to the results of numer- 
ical simulations (e.g., Owocki et al. 1988, Feldmeier et al. 1997b, 
Runacres & Owocki 2002) and observations (e.g., Berghofer et al. 
1997; Miller et al. 2002; Sana et al. 2006). 



1 The equations from the Paper I are denoted as l.x, where x is the equation 
number there. 



2.3 Joint calculation of wind model 

In our previous work we split the calculation of wind model to parts 
below and above the critical point, at which the mass-loss rate of 
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our stationary wind models is determined (Castor, Abbott & Klein 
1975). This effectively means that we neglected the influence of 
the region above the critical point onto the region below this point. 
This was justifiable in our previous models with negligible contri- 
bution of X-ray radiation, because the wind close the critical point 
is optically thick in the UV region and optically thin for other wave- 
lengths (except lines). The splitting of wind solution may not be le- 
gitimate in the case where additional sources of X-rays are present, 
because the X-rays that predominantly originate in the wind above 
the critical point, may penetrate downwards and, consequently, in- 
fluence the solution below this point. 

Here we calculate the joint wind model below and above the 
critical point already during the procedure of mass-loss rate deter- 
mination. To do so, we first calculate several global iteration steps 
(during which we search for the base density corresponding to the 
solution which smoothly passes through the critical point, see Pa- 
per I) for models in which we account only for the wind close 
to the stellar surface. As soon as the appropriate base density is 
roughly known (with a precision of about 30%), we add the solu- 
tion above the critical point and perform several additional global 
iteration steps until the base density is known with sufficiently high 
precision of 1%. This approach enables us to properly take into 
account the influence of X-rays on the mass-loss rate. 



Table 1. Atoms and ions included in the NLTE calculations. Here 'Level' 
means either an individual level or a set of levels merged into a superlevel. 



Ion Levels 


Ion Levels 


Ion Levels 


Ion Levels 


Hi 


9 


Oil 


50 


Alii 


16 


Arm 


25 


HII 


1 


Om 


29 


Aim 


14 


Ariv 


19 


He I 


14 


OIV 


39 


Aliv 


14 


Arv 


16 


Hell 


14 


OV 


14 


AlV 


16 


ArVI 


11 


Hem 


1 


OVI 


20 


AlVI 


1 


ArVII 


1 


ClI 


14 


OVII 


1 


Sill 


12 


Call 


16 


Cm 


23 


Nell 


15 


Sim 


12 


Cam 


14 


Civ 


25 


Nem 


14 


Si IV 


13 


Caiv 


20 


Cv 


11 


Neiv 


12 


Siv 


15 


Cav 


22 


Cvi 


1 


Ne V 


17 


Si VI 


1 


Cavi 


1 


No 


14 


NeVI 


11 


P III 


16 


Fein 


29 


Nm 


32 


NeVII 


1 


P IV 


17 


Feiv 


32 


NIV 


23 


Nail 


13 


P V 


21 


Fe V 


30 


Nv 


13 


Nam 


14 


P VI 


14 


Fe VI 


27 


Nvi 


15 


NalV 


18 


P VII 


1 


Fe VII 


1 


NVII 


1 


Na V 


16 


S II 


14 


Ni III 


36 






NaVI 


1 


S III 


10 


Ni IV 


38 






Mgm 


14 


S IV 


18 


Niv 


48 






MglV 


14 


S V 


14 


Ni VI 


1 






Mg V 


13 


S VI 


16 










Mg VI 


1 


S VII 


1 







2.4 Statistical equilibrium equations 

The statistical equilibrium equation for level i of a given atom (ion) 
has the form (Mihalas 1978, also (1.1)) 



(6) 



where Ni, Nj are relative occupation numbers of studied levels 
(Ni = ni/n atom , where m is the number density of atoms in given 
excitation and ionization state i and n atora is the total number den- 
sity of given atoms; similarly for the level j). Pij are rates of all 
processes by which an atom can change its state (see e.g., Mihalas 
1978). Note that for simplicity we do not explicitly write the radial 
dependence of all variables in equation (6). 



2.4.1 Modification of atomic data 

Model atoms are based on set of TLUSTY files (Hubeny 1988, 
Hubeny & Lanz 1992, Hubeny & Lanz 1995, Lanz & Hubeny 
2003). The original set described in Paper I was extended to allow 
for consistent inclusion of X-rays (see Table 1). Our model atoms 
are based on the data derived from the Opacity Project (Seaton 
1987, Fernley et al. 1987, Luo & Pradhan 1989, Sawey & Berring- 
ton 1992, Seaton et al. 1992, Butler et al. 1993, Nahar & Pradhan 
1993) and Iron Project (Hummer et al. 1993, Bautista 1996, Nahar 
& Pradhan 1996, Zhang 1996, Bautista & Pradhan 1997, Zhang 
& Pradhan 1997, Chen & Pradhan 1999). For phosphorus we em- 
ployed data described by Pauldrach et al. (2001). 



To include the Auger ionization we inserted Auger photoion- 
ization terms into statistical equilibrium equations (6), 

j^i j^i j>i 

where Rfp"' is the Auger photoionization rate with j correspond- 
ing to the ground level of higher ions. Introducing the notation 
ion(i), which means the ionization state, to which the level i be- 
longs, the Auger photoionization rate is given as a product 



(8) 



R Auger _ pAuger 
n ij _ a ion(i)ionO')- rl ion(i)' 

where ai n(i)ion(j) is the Auger yield, i.e., the probability that 
ion(j) — ion(i) electrons are expelled due to the Auger ionization of 
ionic state ion(i), i.e. that ionization state ion(j) is created during 
the process of Auger ionization, and R^j^ is the total inner-shell 
photoionization rate of ion(i). The total inner-shell photoioniza- 
tion rate is given as a sum of partial inner-shell photoionization 
rates from all closed inner shells. Auger rates for transitions ending 
in the ionization states which are not included in statistical equi- 
librium equations are assumed to contribute to the closest lower 
ionization state considered. The term corresponding to the Auger 
ionization is also included in the absorption coefficient in the con- 
tinuum part of the radiative transfer equation. 

The influence of Auger ionization on the temperature and the 
photon emission due to Auger ionization are neglected. 

Photoionization cross sections from individual inner-shells 
were taken from Verner & Yakovlev (1995, see also Verner et al. 
1993) and Auger yields were taken from Kaastra & Mewe (1993). 



2.4.2 Auger ionization 

In addition to the direct ionization (of a valence electron), the Auger 
ionization can significantly alter the ionization state of a stellar 
wind in the presence of X-rays (see Cassinelli & Olson 1979, Olson 
& Castor 1981, MacFarlane et al. 1994, Pauldrach et al. 1994). 



2.4.3 Accelerated lambda iterations 

Iterative solution of the radiative transfer equation together with 
the statistical equilibrium equations may in some cases cause nu- 
merical problems with convergence (see Hubeny 2003 for a re- 
view). To avoid these problems, we included accelerated lambda 
iterations into our models. Their inclusion is based on the method 
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proposed by Rybicki & Hummer (1992). The linearization of 
derived equations is based on Newton-Raphson iterations. Thus, 
our method resembles approximate Newton-Raphson method of 
Hempe & Schonberg (1986). The statistical equilibrium equation 
(6) have for each i the form of 



5>M t£ 
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+ Jv 
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+ J V 



e kT <= dv- 



j<i j 
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-j^- J v du + ^2 N jPji ~ N *J2 P 'i = °' 

(9) 

where we explicitly write rates of direct radiative ionization (I.7a) 
and recombination (I.7b), a;,,, is the photoionization cross-section 
from level i with threshold frequency Vi (similarly for Qj>), J v is 
the mean continuum intensity, T c is the electron temperature, aster- 
isk denotes an LTE value and P„ are rates of all remaining transi- 
tions (collisional ionization, recombination, excitation and deexci- 
tation (1.8) and radiative bound-bound transitions, equations (1.3), 
(1.5), (1.14), and the Auger ionization (8)). 
Introducing the quantity Uji as 

'Ni 



ji — HH^atom [ ~J^~ 



2hi> 6 



hv 
kT,, 



(10) 



the free-bound emissivity has the form of 

V v = N i U H> (") 

where tih is the hydrogen number density and z a tom is abundance 
(number densities ratio) of a given atom relative to the hydrogen. 
During the process of the solution of the radiative transfer equa- 
tion we derive the mean intensity J v from the emissivity (11) as 
follows. Using vector quantities J v = (Jl, Jv, ■ ■ ■ , Jv R ) T and 
j] v — (rjl,r]v, . . . ,r/^ RS ) T , where NR is considered number of 
depth points, the radiative transfer equation for a given frequency 
can be expressed in a symbolic form as 

Jv = *„[rjj, (12) 

where the matrix operator >]/„ represents solution of the radiative 
transfer equation. However, the actual process of solution is differ- 
ent, we do not apply the operator on rj v , but the mean intensity 
J ' v is derived as the solution of linear set of equations 

Vv=^v 1 [Jv], (B) 

which is in fact the formal solution of the radiative transfer equa- 
tion. Hence, during the solution of the radiative transfer equation 
(13) we do not know an explicit form of but we know its in- 
version ^Z 1 ■ Since for the acceleration of convergence of statisti- 
cal equilibrium equations (6) together with the continuum radiative 
transfer equation (13) we need to know derivatives dJv/dNj, we 
consequently need to know the explicit form of >]/„. However, be- 
cause for the solution of equations of statistical equilibrium we use 
only derivatives of mean intensity with respect of Nj at a given 
depth point, we need to know only the diagonal part of the op- 
erator (see equation (11)). This significantly reduces the nec- 
essary computer power. Note that because we use LAPACK pack- 
age (http : / /www . cs . Colorado . edu/ ~ lapack, Anderson 



et al. 1999) for the solution of radiative transfer equation (13), 
which is based on LU decomposition, we can easily calculate di- 
agonal elements of \P„ (see Appendix A). Finally, derivatives of 
mean intensities inserted into the statistical equilibrium equations 
have at the depth point IR the approximate form of 



dJy 
dNi 



m,mUij, 



(14) 



where ^ v ,m,m is the corresponding diagonal element of ^f v at a 
depth point IR and Uij was defined in equation (10). 

During the iterative solution of statistical equilibrium equa- 
tions we solve for the corrections 5N to the actual relative occu- 
pation numbers Ni. These corrections are calculated using equa- 
tion (1.22). Here we also add corrections due to the dependence of 
the mean intensity on Ni (equation 14). Hence, instead of equa- 
tions (1.22) we solve (compare with equation 9) 

Yl^Ni / -^K,m,mUi j SN i du+ 

+ (^) ' £ j£-nm,mU ii 6N j e-*k dv- 

-J2 47vN *{Wj J -^■Km,mU ii SN i e-T* du- 

-^-KjRJRUjiSN,- dv+ 

+E 



N 



dPi* sat , I r> , jit 

31 + Wj dNi 



dNi 



5N 3 



E 



Pij + Ni^-) SN + Ni^-SN 



ON 



dNi 



0, (15) 



where first four sums represent the term of accelerated lambda iter- 
ations and the meaning of other terms is the same as in Paper I. 

Finally, for the calculation of continuum mean intensities J v 
(see equation 13) we use the formal solution (i.e. for given opac- 
ity and emissivity) of the momentum form of the radiative transfer 
equation with inclusion of sphericity factors (Auer 1971). 

To accelerate the convergence even more, we have also in- 
cluded the Ng acceleration (Ng 1974, see also Auer 1987 or 
Hubeny 2003). 



3 STUDIED STARS 

This paper studies the influence of X-rays on the basic parameters 
of radiatively driven stellar winds. However, to avoid using of ad 
hoc stellar parameters, we rather decided to choose a set of pa- 
rameters corresponding to real stars with basic parameters already 
known to some degree of accuracy. 



3.1 Stellar parameters 

Our study is based on O stars with effective temperatures T e g < 
40 000 K, which were detected by ROSAT satellite as X-ray sources 
(Berghofer et al. 1996, hereafter BSC). They were selected out of 
stars studied by Repolust et al. (2004, hereafter R04), Markova et 
al. (2004, hereafter M04), and in the Paper I. To enlarge our dataset 
we also included such stars from Martins et al. (2005, hereafter 
M05) sample, for which measured X-ray fluxes are available in the 
literature. All but one of these stars exhibit the so-called "weak 
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Table 2. Stellar parameters of selected O stars. Stars exhibiting the "weak 
wind problem" appear in the bottom part below the horizontal line. Spectral 
types are taken from the Simbad database. 



oLdT 




Sp. 


D 
it* 


j\/f 

IV 1 


J-eff 


Source 




number 


type 


l K ©J 


l M 0j 


l K l 




^ Per 


24912 


07.5IIIe 


14.0 


36 


35 000 


R04 


a Cam 


30614 


09.5Iae 


27.6 


43 


30 900 


LSL 


AOriA 


36861 


08 III 


12.3 


30 


36 000 


LSL 




54662 


07III 


11.9 


38 


38 600 


M04 




93204 


05V 


11.9 


41 


40 000 


M05 


COph 


149757 


09V 


8.9 


21 


32 000 


R04 


63 Oph 


162978 


08III 


16.0 


40 


37100 


LSL 


68 Cyg 


203064 


08e 


15.7 


38 


34 500 


R04 


19Cep 


209975 


091b 


22.9 


47 


32 000 


R04 


ACep 


210839 


06Iab 


19.6 


51 


38 200 


LSL 


AE Aur 


34078 


09.5Ve 


7.5 


20 


33 000 


M05 


jiCol 


38666 


09.5V 


6.6 


19 


33 000 


M05 




42088 


06.5V 


9.6 


31 


38 000 


M05 




46202 


09V 


8.4 


21 


33 000 


M05 



wind problem", i.e., their theoretically predicted mass-loss rates are 
significantly higher than those derived from observations. 

Adopted parameters of studied O stars are given in Table 2. Ef- 
fective temperatures and radii are taken from R04, M04, M05, and 
Lamers et al. (1995, hereafter LSL). Parameters derived by R04, 
M04, and M05 were obtained using blanketed model atmospheres, 
i.e., they are more reliable than the older ones. Stellar masses were 
obtained using evolutionary tracks either by us (using tracks cal- 
culated by Schaller et al. 1992) or by LSL or M05. The use of the 
evolutionary masses instead of the spectroscopic ones may cause a 
systematic shift due to the well-known discrepancy between these 
masses (e.g. Herrero et al. 1992). However, for many stars from our 
sample these masses are nearly the same. 



3.2 Wind parameters derived from observations 

Wind parameters of studied stars (both derived from observations 
and predicted ones) are given in Table 3. 

X-ray luminosities, which are assumed to originate in the 
wind, were taken from BSC with an exception of stars HD 42088, 
HD 46202, and HD 93204, for which the X-ray luminosities were 
taken from Chlebowski & Garmany (1991) and Evans et al. (2003). 

Because there is still no broad consensus about the influ- 
ence of clumping on mass-loss rates derived from observations, 
we used mass-loss estimations from Puis at al. (2006b, hereafter 
P06), which were regarded as upper limits with respect to clump- 
ing, supplemented by results of M05. Although M05 derived mass- 
loss rates of some stars with inclusion of clumping, for most stars 
selected by us from this sample the clumping factor was set by 
them to one, because these stars exhibit the "weak wind problem". 
The only exception is the star HD 93204, for which M05 provide 
mass-loss rate with clumping taken into the account. However, to 
keep our sample more compact, we used mass-loss rate uncorrected 
for clumping (calculated as M / V/oo, where foo is the clumping 
factor in the outer wind derived by M05). Because P06 concluded 
that mass-loss rates derived from radio data are less influenced by 
clumping, for the star HD 149757 we adopted such rates derived 
by Lamers & Leitherer (1993). 

Terminal wind velocities were taken from LSL, Puis et al. 
(1996), and M04 (the uncertainties were either taken from LSL or 
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Figure 1. The relation between Lx and M plotted using wind parameters 
derived from observations and from theoretical models (fx = 0.02). The 
relations Lx ~ M for stars with optically thick winds in the outer regions 
(stars with large mass-loss rates) and Lx ~ M 2 , valid for stars with op- 
tically thin winds (OC) are also plotted in the graph (dotted lines). Gray 
symbols denote values for stars exhibiting "weak wind problem". 



calculated assuming 10% errors as suggested by Puis et al. 1996). 
For stars that exhibit the "weak wind problem" the terminal veloc- 
ities derived from observations may be just lower limits, similarly 
as in Martins et al. (2004). Consequently, we did not consider the 
terminal velocities of these stars in the following analysis. 



4 THE EMERGENT X-RAYS 
4.1 The L x /L relation 

From the analytic considerations Owocki & Cohen (1999, here- 
after OC) showed that the optically thin X-ray luminosity depends 
on the square of the mass-loss rate Lx ~ (M/vao) 2 , whereas the 
X-ray luminosity of the optically thick wind scales linearly with 
the mass-loss rate Lx ~ M/voo. A slightly different form of these 
relations can be expected due to the assumed dependence of the 
shock temperature Tx on the wind terminal velocity. Consequently, 
in the following we use simpler relations Lx ~ M 2 for the opti- 
cally thin wind in the X-ray region and Lx ~ M for the optically 
thick wind. From the relation M ~ L 1/a ' (Kudritzki & Puis 2000), 
where a' = a — 8 and a and <5 are usual CAK force multipliers, 
the X-ray luminosity of optically thick wind is predicted to be pro- 
portional to L x ~ L 1/a ' (a' « 0.5). 

However, from the observations a slightly different trend 
emerges, because the observed X-ray luminosity Lx is roughly lin- 
early proportional to the total luminosity L (e.g. Chlebowski et al. 
1989; Berghofer et al. 1996; Sana et al. 2006). The origin of the 
difference is not clear. OC suggested that the observed relation can 
be reproduced assuming that fx decreases with radius. 

To understand the origin of Lx/L relation we calculated wind 
models with a fixed filling factor. To be specific, we chose fx = 
0.02. In agreement with the calculations of OC, the X-ray luminos- 
ity of individual stars depends mostly on the wind mass-loss rate 
(see Fig. 1). Stars with large mass-loss rates M > 10~ 7 Mq yr _1 
have an optically thick wind in the X-ray domain (between the out- 
ermost model point and the region where most of the X-rays are 
generated) and the X-ray luminosity of theoretical models calcu- 
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Table 3. Wind and X-ray parameters of studied O stars. Wind parameters derived from observation are described in Sect. 3.2. Predicted wind parameters were 
derived using our NLTE models with different properties of X-ray sources (models with fx = 0, fx = 0.02 and u re i = 0.3). 
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Figure 2. The relation between the X-ray luminosity Lx and the total lu- 
minosity L for studied stars calculated assuming fx = 0.02. Overplotted 
are the average observed relations derived by Berghofer et al. (1997), Sana 
et al. (2006), observed values for studied stars, and the linear fit to the theo- 
retical expectations (only for stars with thick winds). Gray symbols denote 
values for stars exhibiting "weak wind problem". 



lated with a fixed fx depends on M linearly. On the other hand, 
stars with low mass-loss rates M < 10 -7 Mq yr _1 have optically 
thin wind and their theoretical X-ray luminosity depends on M 2 . 

The Lx — M relation plotted using wind parameters derived 
from observations (Fig. 1) is different. The correlation between 
the observed mass-loss rate and X-ray luminosity is less apparent. 
Moreover, the observed X-ray luminosities for stars with low mass- 
loss rates M < 10 _7 Mq yr _1 lie considerably above the theoreti- 
cal expectations. These stars have nearly the same Lx as stars with 
high mass-loss rates. 

The predicted dependence of Lx on L displayed in Fig. 2 re- 
flects the dependence of Lx on M. For stars with dense winds, i.e., 
for stars with large luminosities (L > 5 • 10 38 ergs -1 ) the X-ray 
luminosity scales as M, and due to the dependence of M on L, 
Lx depends mostly on the stellar luminosity as Lx ~ L 1 ' 7 . This 



is in agreement with scaling Lx ~ L 1 ^" where a' « 0.6. For 
stars with thick winds these models roughly recover the observed 
Lx/ L relation. Remaining difference between predicted slope of 
Lx/ L relation and the observed one may be either due the radial 
dependence of the filling factor (as proposed by OC), or due to the 
dependence of the filling factor on wind density (Sect. 8.1). Higher 
scatter of the observed X-ray luminosities may be partly caused by 
poorly known distances (cf. with the results of Sana et al. 2006). 

For stars with lower luminosities, L < 5 • 10 38 ergs -1 , there 
is an abrupt change in Lx/L relation derived from our models with 
fixed fx. These stars have optically thin wind, and their X-ray lu- 
minosities scale as Lx ~ M 2 . Because the mass-loss rate depends 
mainly on the stellar luminosity, the Lx/L relationship derived 
from theoretical models consequently steepens. However, the ob- 
served X-ray luminosities are much higher than the predicted ones. 
The observational Lx/L relation is almost the same for both groups 
of stars with optically thick and thin winds. Because in the case 
of optically thin winds we observe basically all X-rays emitted, 
our description of shock properties of these stars is likely oversim- 
plified. Other possibilities, like different heating mechanism (e.g., 
wind frictional heating, Kubat et al. 2004) are also not ruled out. 

The enhanced X-ray activity may not be limited only to the 
stars exhibiting the "weak wind problem". From observations of 
Berghofer et al. (1997) follows that also many other late O main 
sequence stars show X-ray luminosity corresponding to the stars 
with thick winds. 



4.2 The radial distribution of X-ray emissivity 

Most of X-rays in our models are generated close to the star, as 
shown in the plot of radial variations of X-ray emissivity Fig. 3. 
Due to assumed dependence of the shock temperature on wind ve- 
locity (via equations (1), (3)), the X-ray emissivity is increasing 
function of radius close to the star. On the other hand, we assume 
that the X-ray emissivity depends on the square of the wind den- 
sity (see equation (4)), consequently the emissivity decreases in the 
outer regions. 

For stars with optically thick wind only a very small part of 
generated X-rays finally escapes the wind and may reach a distant 
observer, as can be seen from the difference between the radii at 
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Figure 3. The radial dependence of the the total shock emissivity (integrated 
over frequencies and plotted as 47rr 2 r)x) and the monochromatic luminos- 
ity 4nr 2 F(\) for fx = 0.02 for stars A Cep and £ Per. The quantities are 
plotted relatively to their maximum values. Note that the integrated emis- 
sivity has maximum below the radius at which the X-ray luminosities reach 
the observed value. This indicates that most of X-rays emitted in the wind 
do not escape. 
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Figure 4. The relation between the radius R of the formation of fir lines 
of He-like ions derived from observations by Waldron & Cassinelli (2007) 
and the radius R1/2 at which the monochromatic luminosity is equal to 
half of its terminal value (for fx = 0.02). The radii derived for individual 
lines for giants (III) and main-sequence stars (V) are plotted using different 
symbols. Dashed line denotes one to one relation. 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 



-0.6 
10 





• 'a'6>8 -°° c° * 











. 











f x = 0.02 




BSC • 



10™ 
L [erg s" 



10" 



Figure 5. The comparison between observed (derived by BSC for O stars) 
and model hardness ratio defined by equation (16). 



which 4irr 2 r]x has maximum and at which the monochromatic lu- 
minosities approach the terminal value in Fig. 3. Due to the contin- 
uum X-ray absorption also the X-ray luminosity of optically thick 
wind depends less steeply on the wind mass-loss rate than the X-ray 
luminosity of optically thin wind (see Fig. 1 and OC). The X-ray 
flux is not always monotonically increasing, and in some regions 
it is even negative due to the transfer of part of emitted X-rays to- 
wards the stellar surface. 

In the outer regions the monochromatic luminosities 
47rr 2 F(A) in Fig. 3 rise and reach the value measured by a distant 
observer as the wind becomes optically thin. Waldron & Cassinelli 
(2007) showed that there is a correlation between the radius R^ of 
the formation of forbidden intercombination resonance (fir) lines 
of He-like ions derived from observations, and the radius at which 
the continuum optical depth at associated wavelength is unity. Note, 
however, that the formation radius of the observed X-rays depends 
not only on opacity (optical depth), as Waldron & Cassinelli (2007) 
stressed, but naturally also on emissivity. Consequently, we argue 
that the relationship found by Waldron & Cassinelli (2007) is more 
precisely manifested as a correlation between and the radius 
R 1 / 2 at which the monochromatic luminosity at the associated 
wavelength is equal to half of its terminal value (Fig. 4). 



4.3 The hardness ratio 

To test whether the derived spectral energy distribution of X-rays 
is at least roughly correct, we compare the predicted and observed 
hardness ratio, introduced by BSC as 

HR = (H — S)/(H + S), (16) 

where H and S are the integrated fluxes in the hard (0.5 — 2.0 keV) 
and soft (0.1 — 0.4 keV) energy bands, respectively (see Fig. 5). 
The model hardness ratio is nearly constant for stars with large lu- 
minosities. It decreases for stars with lower luminosities, because 
the stellar wind of these stars is less opaque and regions close to the 
star (where X-rays with lower energies are generated) become vis- 
ible. The observed values for many stars are, however, lower than 
the predicted ones and display larger scatter. In the frame of our 
model this means that the observations cannot be explained just by 
a single value of u rc i. 



5 INFLUENCE OF X-RAYS ON THE IONIZATION STATE 

Although not directly observable, ionization structure of the wind 
gives us a useful direct insight to the effect of different radiation 
processes in the stellar atmosphere and wind. In our models the X- 
rays affect wind ionization structure via direct and Auger ionization 
(MacFarlane et al. 1994, Pauldrach et al. 1994). 
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Figure 6. Radial variation of nitrogen ionization fractions in the wind model 
of A Cep for fx = 0.0028 that fits the observed Lx- Top panel: Compar- 
ison of models calculated with additional X-ray sources (solid lines) and 
without these sources (dashed lines). Bottom panel: Comparison of models 
with additional X-ray sources calculated with inclusion of Auger processes 
(solid lines) and neglecting the Auger processes (dashed lines). For better 
readability, the ionization fraction of N II is not plotted. 



5.1 Radial variations of the ionization fractions 

As an example we discuss the variations of nitrogen ionization in 
the case of the star A Cep (HD 210839). To understand the influence 
of individual ionization mechanisms on the wind ionization struc- 
ture we also calculated models without additional X-ray sources 
and models which take these sources into account, but which ne- 
glect Auger processes (see Figs. 6 for the plot of the variations of 
nitrogen ionization fractions with the radius). 

The ionization fractions of most abundant ions (NlII, NlV) 
close to the stellar surface are not significantly influenced by the 
additional sources of X-rays. This region is optically thick for the 
X-ray radiation and only a very tiny amount of X-rays emitted in 
the highly supersonic wind parts may penetrate close to the stel- 
lar surface and modify the ionization fractions of minor higher 
ions there (in our case N VI and N VII). With increasing radius (for 
r/R f — 1 > 0.1) the wind becomes more ionised due to lower 
wind density, and X-rays start to play a more important role in the 
ionization balance. As a result of our assumptions, the Auger ion- 
ization significantly influences the ionization balance only in the 
outer wind regions. This is caused by the fact that more energetic 
X-rays emitted in the wind with sufficiently high velocity are nec- 
essary for the Auger ionization than for direct ionization. 

The ionization fraction of N V is increased due to direct ion- 
ization from N IV and is influenced also by ionization to and recom- 
bination from N VI. The ionization fraction of N VI is increased by 




Figure 7. The radial dependence of selected sums of radiative ioniza- 
tion/recombination rates NiPij in the model of A Cep (fx = 0.0028). 
Rates are summed over all transitions between corresponding ions. Close 
to the star, the ionization fraction of N V is given by the radiative ionization 
balance with N IV. In the outer regions also the recombination from N VI in- 
fluences the fraction of N V . Consequently, N V partly occurs due a two-step 
process (Auger ionization of N IV followed by the radiative recombination). 



both Auger and direct ionization. N VII is generated due to the di- 
rect ionization of N VI (see Fig. 7). 

The result that for r/R* — 1 < 0.5 the ionization fraction of 
N VII is higher in the case of neglected Auger processes than in the 
model with these processes may seem paradoxical, but it has a sim- 
ple explanation. In the case of neglected Auger processes the stellar 
wind is less opaque in the X-ray region and, consequently, the in- 
tensity of X-ray radiation is higher. Because N VII is created via 
direct ionization of N VI, the ionization fraction of N VII is higher. 



5.2 Variations with the effective temperature 

To compare the ionization equilibrium for different stars, we plot 
the ionization fractions of selected ions for the points where the 
wind velocity is equal to v r = 0.5«oo (Fig. 8). The model predic- 
tions are compared with the ionization fractions derived from ob- 
servations by Howarth & Prinja (1989), Lamers et al. (1999) for 
Galactic stars, and Massa et al. (2003) for stars from the Large 
Magellanic Cloud. Stars from the Clouds have generally lower 
metallicity than the Galactic ones (e.g. Martins et al. 2004). How- 
ever, because the dependence of ionization fractions on metallicity 
is not strong, we use also these fractions for our comparison. 

The ionization fractions are influenced not only by the stellar 
effective temperature and by local wind parameters, but also by the 
amount of additional ionizing X-rays. Consequently, the variations 
plotted in Fig. 8 are far from being monotonic. In agreement with 
MacFarlane et al. (1994), the ionization fractions of the dominant 
ions are usually not influenced by X-ray emission. X-rays above 
all modify the ionization fractions of minor higher ionization states 
of a particular atom both due to the direct and Auger ionization. 
There is a satisfactory agreement between our results and models 
of MacFarlane et al. (1994). Note the influence of X-rays on the 
ionization balance already very deep in the wind as both observa- 
tional and theoretical results show relatively large amount of O VI 
close to the stellar surface for v r = 0.05 Woo (Fig. 9). However, 
X-rays do not alter the ionization state of the dominant and of the 
lower neighbouring ions there (cf., Fig. 6). 

The influence of X-rays on the ionization structure is stronger 
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Figure 8. Ionization fractions as a function of the effective temperature for individual stars from our sample at the point where the radial velocity v r = 0.5i>oo • 
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for stars with less dense winds and for stars with weaker flux of the 
ionizing radiation, i.e. for cooler stars. 

The presence of N V lines in the wind spectra is a frequently 
used argument for the existence of an additional source of ioniza- 
tion. However, as we can see from the comparison of the ioniza- 
tion structure derived from observations and from our models for 
v = 0. 05«oo, and v = 0.5i>oo (Figs. 9, 8), the ionization fraction of 
N V there can be explained just by the detailed NLTE models with- 
out a contribution of X-ray radiation. This was shown already by 
Pauldrach (1987). Only the presence of N V in a significant amount 
at higher wind speeds (see Fig. 10 and also Fig. 6) is caused by the 
existence of a strong X-ray source. 



5.3 Phosphorus 

The failure of present NLTE wind models to reproduce the ion- 
ization fractions of Pv derived from observations is one of the 
main arguments for lower values of wind mass-loss rates. There 
is large discrepancy between P V ionization fraction derived by us 
and those derived from observations by Fullerton et al. (2006), who 
concluded that either on average the ionization fraction of P V is 
lower than 0.1 or there is a significant disagreement between mass- 
loss rates derived from P V and Ha lines. Some part of this discrep- 



ancy may come from the neglected influence of inhomogeneities on 
the wind ionization fractions or radiative transfer (Puis at al. 2006a; 
Krticka et al. 2007; Oskinova et al. 2007). Anyway, it is worthwhile 
to understand whether some part of this discrepancy is not caused 
by the fact that P V is a fragile ion, which is not a dominant one for 
any effective temperature, and to test the influence of X-rays on P V 
ionization fraction. 

Our calculations show that Pv is a dominant ion in the in- 
ner part of the wind only for stars with T e ff < 34 000 K (in accor- 
dance with theoretical results discussed by Fullerton et al. 2006, see 
right lower panels of Fig. 8 for the the predicted ionization fraction 
of P V and Piv). For hotter stars also P VI becomes an important 
ionization stage. However, the presence of shock X-ray radiation 
does not significantly modify the ionization fraction of P V. Conse- 
quently, even for hotter stars for which P V is not a dominant ion 
there is a significant discrepancy between theoretical and observa- 
tional results, as can be seen from Fig. 11. Here we plot the values, 
which can be in fact derived from observations, i.e. the product of 
the ionization fraction and the mass-loss rate qp v M, and compare 
it with the predicted ones. For hotter stars (T e ff > 35 000 K), the 
predicted values are on average by a factor of 7 higher than those 
derived from observations (with an exception of the star HD 42088, 
which displays the "weak wind problem"). On the other hand, for 
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Figure 12. The influence of X-rays on the radial variations of phospho- 
rus ionization fractions in the wind model of the star HD 203064. Note 
that X-rays do not significantly change P V ionization fraction and that P V 
dominates for r > 2R* . 



1.2 



10" 



10" 



10" 



10"' 



10"' L 



• MIM(L=0) 
' v„ I v„{f x =0) 



. .. 



10"' 
10" : 
10" : 
10"' 
10" E 
10" £ 



•+ 


• 
• 


• 

• 


-0 








o 


+ o 




+ 







o 


o 


+ 










+ 














+ 


o o ° 










o 






o 








o 












o 
















N V 



0.9 
10" 



10" 



30000 



35000 



40000 



Figure 10. The same as Fig. 8, however for v r = 0.8i>oo- 



10 J 



I 1 



•5 10 1 

> 



P V 



1 ' — ■ 

30000 



35000 



40000 



Figure 11. Comparison of the product of the P V ionization fraction and the 
wind mass-loss rate derived from observations (Fullerton et al. 2006) and 
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the star HD 42088, that exhibits the "weak wind problem". 
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Figure 13. Variations of the mass-loss rate M and the terminal velocity 
Voo of a Cam with fx and Lx (£x depends linearly on fx) relative to their 
values derived using the models with no shock X-ray sources. 



cooler stars (T e ff < 35 000 K) the predicted values are on average 
even about 100 times higher than those derived from observations. 

The maximum frequency considered in our models (corre- 
sponding to 1 .7 keV) is lower than P V inner-shell ionization energy 
threshold 2.2 keV. To test our results, we calculated also models 
with higher maximum frequency considered, i.e., corresponding to 
3.7 keV (see Fig. 12). However, even in these models the Auger 
ionization of P V due to low flux of high energy X-rays does not 
significantly modify the ionization fraction of P V. 



6 INFLUENCE OF X-RAYS ON HYDRODYNAMIC 
STRUCTURE 

To understand the influence of X-rays on wind parameters, we cal- 
culated several wind models corresponding to a Cam stellar pa- 
rameters, however with a different value of fx, i.e. with a different 
value of X-ray luminosity. The variations of wind parameters with 
X-ray luminosity are shown in Fig. 13. As a consequence of as- 
sumed dependence of X-ray emissivity on fx via equation (4), the 
emergent X-ray luminosity depends linearly on fx- For lower X- 
ray luminosities (Lx/L < 10~ 8 ) both mass-loss rate and terminal 
velocity are not significantly affected by the X-ray radiation, since 
the latter influences only the ionization state of minor ions (typi- 
cally VI-VII). This picture changes for higher X-ray luminosities 
(Lx/L > 10~ 8 ). The mass-loss rate slightly decreases with in- 
creasing X-ray luminosity, whereas the terminal velocity slightly 
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Figure 15. X-ray continuum opacity calculated with and without shock X- 
ray emission for £ Per. Here we give the opacity per unit of mass averaged 
for the radii 1.5Rq < r < 5 Rq . The inclusion of X-rays leads to a slight 
decrease of the opacity in the X-ray region. 



increases. Changes of wind parameters are caused by modification 
of ionization equilibrium due to X-rays. The X-ray radiation shifts 
the ionization balance to higher excited states of all elements. Be- 
cause the number density of minor ionization state Fe IV, which is 
important for the radiative acceleration close to the star, is lowered, 
the radiative acceleration decreases. Consequently, the wind mass- 
loss rate slightly decreases in the case of a Cam. 

Generally, the changes of the wind mass-loss rate due to X- 
rays are relatively small. The mass-loss rates for the same star cal- 
culated with fx from to 0.04 differ by no more than about 5%. 
This is connected with the fact that in our models the wind mass- 
loss rate is determined in the region close to the stellar surface be- 
low the critical point. Because this wind region is opaque to X-rays, 
the X-rays do not influence the ionization states significantly con- 
tributing to the radiative force there and the mass-loss rates remain 
basically unaltered. 

The change in ionization due to the presence of X-rays does 
not decrease the line-force due to the lighter elements. The impor- 
tant lines stay optically thick and higher excited elements capable 
to contribute to the radiative force emerge (e.g., Nv, Ov, Ovi). 
Thus, the situation in the outer wind, where the contribution of iron 
lines is smaller and the influence of X-rays is higher, is different. 
The increase of the radiative force (supported in some cases by the 
decrease of M) causes the increase of the wind terminal velocity. 

For stars with stronger ionizing flux in the far UV region the 
influence of X-rays on the wind structure is not so significant. More 
X-rays are necessary to modify the ionization fractions of wind 
driving ions. Consequently, the increase of the wind terminal ve- 
locity is larger for cooler stars (T e g < 35 000 K, see Fig. 14) or for 
stars with more tenuous winds. 

The energy of artificially included X-rays also transforms to 
the thermal energy due to the X-ray absorption. Consequently, wind 
temperature may be by few thousands Kelvin higher due to X-rays. 
However, this effect does not significantly influence the radiative 
force. 



7 CONSEQUENCES FOR THE X-RAY LINE TRANSFER 

Using modern X-ray spectrographs it became possible to resolve 
the X-ray line profiles of hot stars. The study of fir lines of He-like 
ions enabled to derive radii at which the X-ray emission is gener- 
ated (e.g., Kahn et al. 2001; Miller et al. 2002; Raassen et al. 2005; 
Leutenegger et al. 2006). It is also possible to predict the shapes of 



these X-ray lines (e.g., Owocki & Cohen 2001). Most of the predic- 
tions assume that the X-ray lines originate in the hot optically thin 
environment surrounded by cool wind, which is optically thick in 
continuum. These theoretical calculations based on currently avail- 
able theoretical mass-loss rates predict asymmetric X-ray line pro- 
files due to the wind continuum absorption. However, the observed 
X-ray line profiles are in most cases symmetric (e.g., Miller et al. 
2002; Leutenegger et al. 2007). This either indicates that the the- 
oretical wind mass-loss rates are overestimated or this leads to the 
discussion on the influence of wind inhomogeneities on the X-ray 
line transfer (Feldmeier et al. 2003; Owocki et al. 2004; Oskinova 
et al. 2006a; Owocki & Cohen 2006). Here we study the influence 
of modified ionization equilibrium due to hot shock emission on 
the line transfer in the X-ray region. 



7.1 Influence of shock emission on the X-ray opacity 

The theoretical prediction of the asymmetric shape of the X-ray line 
profiles is sensitive to the value of the opacity in the X-ray region. 
Krticka & Kubat (2007) concluded that a more realistic (lower) 
value of adopted metallicity leads to lower continuum opacity in 
the X-ray region, and consequently to a better agreement between 
theoretical and observed X-ray line profiles. Here we test whether 
the ionization shift due to the presence of X-rays may affect the 
X-ray continuum opacity. 

Our calculated opacities are in a good agreement with those 
presented by Oskinova et al. (2006a). Inclusion of shock X-ray 
emission leads to a slight decrease (on average by 8%) of the cal- 
culated opacity per unit of mass due to the ionization shift of ions 
He II, C V, N IV, O IV etc., as shown in Fig. 15. This could slightly 
improve the agreement between predicted and observed X-ray line 
shapes, but the decrease is too small to explain the main part of the 
discrepancy. 



7.2 Influence of the resonance-line scattering in the cool wind 
on the line transfer in the X-ray region 

Ignace & Gayley (2002) proposed that X-ray emitting regions may 
be optically thick in lines. This effect can lead to more symmetric 
X-ray line profiles, which better correspond to the observed ones 
(see also Leutenegger et al. 2007). However, the ambient cool wind 
is usually assumed to absorb the X-ray radiation only in continuum 
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Figure 16. The Sobolev optical depth (17) due to the resonance line scatter- 
ing due to Is 2 X S - ls2p l P lines of helium-like ions N VI (A = 28.787 A) 
and O VII (A = 21.602 A) in the cool wind of HD 93204. The model was 
calculated for fx = 0.0094 that fits the observed X-ray luminosity. 



and the line absorption within the cool wind is being neglected. 
Here we test whether this assumption is adequate. 

He-like ions N VI and O VII, whose lines are observed in the 
X-ray spectra of hot stars, are present in a nonneglible amount also 
in a cool ambient wind (see Fig. 8). This is caused by the direct and 
Auger ionization of corresponding less ionised ions. The influence 
of the resonance-line scattering in the cool ambient wind on the 
profiles of X-ray lines emitted in shocks can be estimated from the 
optical depth of these lines in the cool wind. Neglecting the occu- 
pation number of upper level n;/<7» » fij/gj, the Sobolev optical 
depth in a given line is within the radial streaming approximation 
roughly given by (e.g., Castor 1974) 



ne { dv r 

Uijij I — 
meVii V ar 



(17) 



where Vij is the line frequency, fij is the oscillator strength, and 
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Figure 17. The X-ray luminosity calculated assuming the dependence of 
fx on the cooling time via equation (18) in comparison with observations. 
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Figure 18. The radial variations of ex (see equation (18)) for selected stars. 



weights gi, gj. From this equation it follows that the effect of the 
resonance line scattering in the cool wind is effective mainly for 
stars with high mass-loss rates or for stars where X-rays more sig- 
nificantly influence the ionization equilibrium (i.e., stars with large 
fx, where higher fractions of He-like ions emerge). 

For stars HD 93204, 162978, 209975, and 210839 the reso- 
nance line scattering in the cool wind is able to influence the emer- 
gent X-ray line profiles of N VI and in some cases also of O VII (see 
Fig. 16). This effect, discussed also by Oskinova et al. (2006b), is 
sensitive to the properties of X-ray sources. 



8 DISCUSSION 

8.1 The dependence of fx on the cooling time 

The parameter fx is an important free parameter that determines 
the wind X-ray luminosity. As it was used for some models, it is 
possible to specify fx individually for each star using observed X- 
ray luminosity. As with fixed fx we predict steeper Lx — L relation 
than observed (Fig. 2), fx which is inversely proportional to the 
luminosity would lead to a Lx — L relation that better corresponds 
to the observed one. 

The inverse proportionality of fx on L (or on the wind density) 
may be a more realistic one. The total X-ray emissivity is given by 
the amount of energy dissipated by the shocks and the fraction of X- 
ray emitting material (or fx) may be determined by this dissipated 
energy and by the cooling length. In fact, fx may be lower for stars 
with dense winds in which the hot material is able to cool down 



more efficiently. To correct our results for the dependence of fx on 
the cooling time, we calculated wind models where fx is scaled by 
the value of the cooling time relatively to the hydrodynamical time- 
scale (cf. Antokhin et al. 2004; Zhekov & Palla 2007), i.e. roughly 
given by 



fx = /oex = fo 



Oh/3 



(£p) 2 Jo°° A ( Tx ) Av r 



(18) 



where we assume fo = 0.7. Since according to (18) fx ~ r, Lx 
now increases ad infinitum. To avoid a possible divergence of Lx 
we limit fx in (18) by a certain value, namely fx < 0.1 as in the 
outer regions shocks may become adiabatic or weaker (Feldmeier 
et al. 1997a,b). The derived relation between Lx and L in Fig. 17 
shows a good agreement with observations indicating that the vari- 
able cooling time may indeed be an important property of wind 
shocks. 



8.2 The weak wind problem 

Comparison of mass-loss rates derived from observation with the 
predicted ones (Table 3) shows that there is a good agreement be- 
tween observation and theory for stars with large mass-loss rates 
M > 10~ 7 M Q yr _1 (provided that the observations are not signif- 
icantly influenced by clumping). The situation is markedly different 
for stars with lower mass-loss rates M < 10~ 7 M yr _1 (except 
HD 149757). The predicted mass-loss rates are generally by order 
of magnitude higher than those derived from observations. These 
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in the post-shock region is longer. Consequently, fx may become 
higher for lower metallicities. In such a case we can expect the 
dependence of fx on the metallicity roughly as fx ~ Z^ 1 . Within 
this approximation the total X-ray luminosity would scale as 
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Figure 19. The ratio between the calculated X-ray fluxes from models with 
metallicities Z = 0.9 and Z = 1 (for the same fx = 0.02). The dashed 
lines denote the approximate relations (equations (19), (20)). 

stars exhibit the so-called "weak wind problem", whose solution 
is not yet known (Bouret et al. 2003; Martins et al. 2004, M05, 
Krticka 2006). 

The X-ray activity of these stars is higher than that of stars 
which do not exhibit the weak wind problem. Within our models 
this means that their X-ray luminosity predicted assuming a fixed 
filling factor is by two to three orders of magnitude lower than the 
observed one (see Fig. 2). However, the X-ray emission of these 
stars may still originate in the wind, because their X-ray luminosity 
is lower than the wind kinetic energy lost per unit of time l/2Mv% 
(M05). Note that due to low wind density of these stars the cooling 
length may become comparable with the hydrodynamical scale (as 
indicates the plot of ex in Fig. 18). Consequently, there is a possi- 
bility that once the wind of these stars is heated by the shocks it is 
not able to cool down sufficiently, and remains hot and unaffected 
by the radiative acceleration. This can be a reasonable explanation 
of the "weak wind problem" and high X-ray luminosity Lx of these 
stars. A similar scenario was invoked by Cohen et al. (2008) to ex- 
plain the X-ray spectrum of (3 Cru. 

8.3 The metallicity dependence 

To discuss the influence of metallicity on Lx we recalculated our 
models with lower value of the mass fraction of heavier elements 
Z. Because the dependence of the free parameter fx on metallicity 
is not known, we calculated these models with fixed fx- 

The derived variations of Lx with Z given in Fig. 19 can 
be easily interpreted. The mass-loss rate depends on metallicity 
roughly as M ~ Z ' 67 (Krticka 2006). The X-ray emission origi- 
nates mainly from lines of heavier elements, thus we can assume 
that the emissivity per unit of mass depends on the metallicity 
roughly as A„(Tx) ~ Z. Consequently, from the relations of OC 
(see Section 4.1) the optically thin X-ray luminosity scales with the 
metallicity as 

Lx ~ Z 2 ' 34 . (19) 

In the optically thick case we can assume that for higher fre- 
quencies {y > 1 x 10 17 Hz) the X-ray opacity Kx is mostly due 
to direct and inner-shell ionization of heavier elements. As a result 
of this Kx ~ Z, and in the optically thick case the relations of OC 
yield Lx ~ MA„/kx, i.e. 

L x ~ Z - 67 . (20) 

The relations (19) and (20) are derived assuming that fx does 
not depend on the metallicity. However, lower value of A for lower 
metallicities could effectively mean that the cooling time of the gas 



8.4 Multicomponent effects 

The multicomponent effects due to inefficient collisions between 
individual wind particle components may significantly influence 
structure of low-density winds (e.g., Springmann & Pauldrach 
1992; Krticka 2006). The Coulomb frictional force between wind 
components depends on the charge of these components, which 
may be influenced by the presence of X-rays. However, our tests 
showed that such influence is in most cases insignificant, because 
the X-rays typically modify the fraction of minor ionization stages 
only. 

We tested whether the multicomponent effects are important 
for stars exhibiting the "weak-wind problem". For this purpose we 
treated heavier elements as individual components (as in Krticka 
2006). We concluded that phosphorus due to its low abundance may 
decouple from the wind of HD 34078 and HD 38666, and the ve- 
locity difference between phosphorus and hydrogen is high also in 
the wind of HD 46202. For these stars the multicomponent effects 
and the ionization shift induced by it may be important. 

8.5 Model simplifications 

The inclusion of additional X-ray sources into our NLTE wind 
models was done assuming a simplified shock picture, in reality 
these sources may have slightly different properties. For example, 
the X-ray source can not be described just by one temperature 
(Feldmeier et al. 1997a). However, the main result of this paper, 
i.e., that the presence of X-rays does not significantly influence the 
wind parameters of O stars was justified for a broad range of shock 
parameters. On the other hand, the predicted ionization structure is 
much more sensitive to properties of X-ray sources. 

There is a growing observational evidence that hot stars winds 
are clumped (Bouret et al. 2003; Martins et al. 2004; Puis at al. 
2006b). We postpone the effect of clumping on the wind structure 
for a special study. The clumping may influence not only the com- 
parison of ionization fractions derived from observation and from 
the theoretical models, but also modify the radiative force and the 
mass-loss rates (Krticka et al. 2007; de Koter et al. 2007). On the 
other hand, the effect of porosity (Owocki et al. 2004; Oskinova et 
al. 2006a) may lead to decrease of the effective opacity in the X-ray 
region. This may affect the total X-ray luminosity, but as the poros- 
ity affects also the predicted wind parameters (Krticka et al. 2007) 
and as its existence in the wind is a matter of discussion nowa- 
days, we do not include it into our models. Note that the present 
approach of the inclusion of both clumping and porosity into wind 
models is based on free parameters whose values are not very well 
observationally nor theoretically constrained. 



9 CONCLUSIONS 

We studied an effect of additional X-ray emission on our NLTE 
wind models. To this end, we assumed that some part of the wind 
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material is heated to a very high temperature (of order of 10 6 K) 
and emits X-rays, and we studied their influence on wind structure. 

Using our wind models we derived the same scaling of the 
total X-ray luminosity i x with the mass-loss rate M as Owocki & 
Cohen (1999). For stars with optically thick winds most of X-rays 
are absorbed in the wind, and the X-ray luminosity scales with wind 
mass-loss rate roughly as L x ~ M. For stars with optically thin 
wind in the X-ray region the X-ray luminosity scales as L x ~ M 2 . 

The models with fixed filling factor fx can roughly explain the 
observed L x /L relation of stars with optically thick wind in the X- 
ray region. The remaining deviation between the observed slope of 
the Lx/L relation and that derived using fixed fx can be the result 
of the dependence of the shock cooling time on wind density. The 
stars with optically thin wind exhibiting the "weak wind problem" 
emit more X-rays than predicted by a simple optically thin scaling. 
Their enhanced X-ray emission and the "weak wind problem" itself 
may be caused by large wind shock cooling time, comparable with 
the hydrodynamical time scale. 

The X-rays influence the ionization state in the whole wind. 
However, a significant amount of X-rays is emitted only in highly 
supersonic wind regions. These X-rays mostly do not penetrate 
close to the stellar surface, consequently the ionization changes in- 
duced by X-rays there are not large enough to change the radiative 
force there significantly. Because the mass-loss rate is determined 
in the region close to the stellar surface, the presence of X-rays does 
not significantly change its value. On the other hand, the terminal 
velocity is determined in the outer wind regions where the influence 
of X-rays is more important. Consequently, the wind terminal ve- 
locity may be slightly affected by X-rays, especially for stars with 
weaker ionizing continua (with T e tf < 35 000 K). 

Although the Auger processes are important for the ionization 
balance in the outer wind regions, the non-negligible presence of 
higher ionization stages (sometimes also called "superions") is a 
consequence of also a direct ionization. 

We compared the model ionization fractions with those de- 
rived from observations. For higher ionization stages the inclusion 
of X-rays leads to better agreement between theory and observa- 
tions. However, for many ions a significant discrepancy between 
theory and observations still remains. The discrepancy between the 
mass-loss rates derived from the Ha line or radio emission and 
P V line profiles led Fullerton et al. (2006) to conclusion that real 
wind mass-loss rates of hot stars are much lower than the "classi- 
cal" ones. This conclusion is sensitive to the assumed P V ionisation 
fraction. We have shown that for stars with T eff > 34 000 K the P VI 
ion becomes a dominant one, which gives a better agreement for 
hotter stars. But even for these hotter stars a significant disagree- 
ment between theory and observations remains. We did not find 
any significant influence of X-rays on P V ionization fraction. Our 
calculations suggest that the effect of X-rays cannot be the reason 
the P V mass-loss rates are low. This problem is especially severe 
for stars with T eff < 34 000 K for which we predict P V to be the 
dominant ionization stage in the wind. 

We studied the influence of modified ionization equilibrium 
due to shock emission on the X-ray line transfer. While continuum 
absorption by the cool wind in the X-ray domain is not significantly 
affected by X-rays itself, additional line opacity sources in a cool 
wind emerge (see also Oskinova et al. 2006b). This may affect the 
line profiles of more abundant elements (C, N, and O) observed in 
the X-ray spectra. 

We can conclude that although the inclusion of X-rays into 
NLTE wind models leads to a better agreement between theory and 
observation, it can not explain the most painful problems of present 



hot star wind research, i.e., lower mass-loss rates derived from ob- 
servation with taking into account clumping, and too-low ionization 
fractions of P V derived from observations. 
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APPENDIX A: CALCULATION OF DIAGONAL 
ELEMENTS OF 

The equation of the radiative transfer (13) is solved using the LU 
decomposition, 

= LU, (Al) 

where L is a product of permutation and lower triangular matrices 
with 1 subdiagonal, and U is upper triangular matrix with 2 su- 
perdiagonals. In such a case the inversion matrix >]/, y can be easily 
calculated (Ralston 1965) as 

*„ = U _1 L- 1 , (A2) 

where IP 1 and L _1 are upper and lower diagonal matrices. If we 
denote u = IT 1 and 

Ui,< =U^\ (A3) 

_ Uj,j-iUj-i,j + Ui,j- 2 Uj-2,j .... 
U »J — Tj ) l/"7 

U 3 ,3 

and similarly for L _1 . From this, diagonal elements of ^ v can be 
easily calculated. 

Instead of using matrix inversion it is possible to derivate 
(13) with respect to the occupation number Ni and calculate par- 
tial derivatives dJ^/dNj as a solution of derived system of linear 
equations. Although this would be likely fast solution (remember 
that since we employ LU decomposition for the solution of equa- 
tion 13, we can use this solution procedure to solve other equations 
with the same matrix), we would have to solve these equations for 
each Nj, and this would be very time consuming. 
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